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ABSTRACT 


OBJECTIVE  ANALYSIS  OF  A  TWO-DIMENSIONAL  DATA  FIELD 
BY  THE  CUBIC  SPLINE  TECHNIQUE 

A  procedure  to  use  the  spline  interpolation  technique  on  an 
arbitrarily  prescribed  two-dimensional  data  field  is  described.  In 
order  to  use  this  technique  it  is  necessary  to  obtain  an  initial 
approximation  to  the  data  at  the  grid  points .  This  is  achieved  by 
fitting  spherical  surfaces  to  the  data.  Bi-directional  spline  inter¬ 
polation  is  then  applied  repeatedly  on  the  grid  point  estimates  of  the 
data  to  produce  convergence  to  the  true  surface. 

The  spline  interpolation  technique  and  another  objective  analysis 
technique  developed  by  Gilchrist  and  Cressman  are  tested  against  an 
exact  solution  and  the  resulting  analyses  are  compared.  Real  tempera¬ 
ture,  geopotential  height,  and  wind  data  for  various  pressure  surfaces 
are  analyzed  by  the  spline  method  and  the  results  are  compared  to 
subjective  analyses  of  the  same  data. 
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1.  Introduction 


With  the  development  and  continuing  improvement  of  the  electronic 
computer  came  the  beginning  of  practical  numerical  weather  prediction. 

The  present  day  computers  make  it  possible  to  use  elaborate  numerical 
models  to  improve  forecast  accuracy;  however,  there  are  certain  limi¬ 
tations  to  this  improvement  which  cannot  be  overcome  by  increasing  the 
speed  and  capacity  of  the  computer  or  by  improving  the  set  of  equations 
which  define  the  model.  The  two  major  limitations  are:  1)  the  density 
of  the  network  of  observation  stations;  and  2)  the  una\'oidable  error, 
instrument  or  otherwise,  which  is  introduced  while  making  any  type  of 
measurement.  The  first  limitation  imposes  a  well  defined  limit  on  the 
scale  of  atmospheric  motion  which  can  be  resolved  by  the  model.  Of 
course  the  solution  to  this  problem  is  simply  to  increase  the  density 
of  the  observing  network  thereby  obtaining  the  "scale  of  data"  necessary 
for  predicting  the  corresponding  scale  of  motion. 

The  second  limitation  is  far  more  difficult  to  eliminate.  Any  set 
of  observations  contains  certain  random  errors  as  well  as  some  small 
scale  fluctuations  (noise) .  In  using  this  data  in  any  numerical  model 
it  is  of  utmost  importance  to  eliminate  these  errors  and  the  noise  if 
the  final  results  are  to  be  meaningful.  This  has  been  demonstrated,  for 
example,  by  the  studies  of  Best  (1956)  and  Barggran  (1957).  Since 
meteorological  prediction  is  an  initial  value  problem,  it  becomes  impera¬ 
tive  to  begin  any  numerical  weather  prediction  with  the  "beat  possible” 
representation  of  the  real  data. 

Various  attempts  at  solving  this  problem  have  been  introduced  in 
meteorology  within  the  past  fifteen  years  and  shall  be  referred  to  as 
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"objective  analysis".  An  appropriate  definition  of  objective  analysis 
was  given  by  Gandin  (1963). 

"Objective  analysis  includes  the  development  and  realization 
of  methods  which  make  it  possible  to  use  the  measurement 
data  of  meteorological  stations  to  reconstruct  objectively 
the  fields  of  the  meteorological  elements  [variable],  or  at 
any  rate  to  specify  their  values  at  the  nodes  [grid  points] 
of  some  type  of  regular  network. 

Actually,  objective  analysis  includes  three  distinct  functions: 

1)  elimination  or  correction  of  gross  errors  in  the  data  field;  2) 
interpolation  of  data  to  obtain  values  on  a  grid;  and  3)  smoothing  of 
the  resulting  values  at  the  grid  points. 

Probably  the  first  attempt  at  objective  analysis  of  meteorological 
data  was  by  Kibel  in  1949.  Kibel  used  formulas  based  on  the  method  of 
least  squares  to  describe  the  field  of  data  by  second  and  third  order 
polynomials.  Also  in  .1949,  Panofsky  represented  a  field  by  cubic 
polynomials  and  showed  that  the  introduction  of  random  observational 
errors  into  the  data  field  resulted  in  only  minor  variations  in  the 
polynomial  representation. 

By  1954,  the  need  for  a  better  objective  analysis  technique  to  be 
used  in  conjunction  with  the  rapidly  developing  field  of  numerical 
weather  prediction  resulted  in  the  method  developed  by  Gilchrist  and 
Cragsman  (1954).  Their  method  was  based  on  fitting  a  second  degree 
polynomial  by  the  method  of  least  squares  to  the  data  in  a  limited  araa 
around  each  grid  point.  Wind  values  were  incorporated  in  the  scheme 
by  using  the  geos  trophic  assumption  to  determine  *h  (where  h  is  the 
deviation  ol  height  from  the  standard  atmosphere)  at  a  data  point. 
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Thus  each  data  point  supplied  three  pieces  of  information,  h,  and 

'  ’  X 

9h 

to  be  used  in  the  least  squares  fitting.  For  regions  of  little 
data,  this  method  proved  to  be  inadequate  since  at  least  six  initial 
pieces  of  information  were  needed  to  determine  the  second  degree  poly¬ 
nomials  which  defined  the  field.  In  fact,  Gilchrist  and  Cressman  found 
that  with  less  than  ten  pieces  of  data  the  calculation  was  subject  to 
significant  error.  In  regions  of  sufficient  data,  however,  they  found 
that  numerical  predictions  based  on  the  objective  analyses  were  an 
improvement  over  those  based  on  subjective  analyses. 

About  a  year  after  Gilchrist  and  Cressman  (1954)  introduced  their 
objective  analysis  technique,  Bergthorssen  and  Doos  (1955)  developed  a 
new  approach  to  objective  analysis.  Their  method  differed  from 
Gilchrist  and  Cressman's  in  that  they  first  determined  what  is  now 
called  the  "preliminary  field".  This  field  was  computed  from  the 
weighted  mean  of  the  forecast  values  at  the  grid  points  and  the  clima¬ 
tological  norms  for  that  time.  The  data  obtained  from  the  observation 
stations  was  then  used  in  a  sequence  of  three  corrections  which  were 
then  applied  to  the  preliminary  field  with  different  weights,  depending 
on  the  distance  from  station  to  grid  point.  The  numerical  predictions 
based  on  the  objective  analyses  technique  of  Bergthorssen  and  Doos  were 
approximately  the  same  as  the  predictions  based  on  the  subjective 
analyses , 

Approximately  five  years  after  the  development  of  his  first  objec¬ 
tive  analysis  technique,  Cressman  (1959)  introduced  a  modification  of 
the  Bergthorssen-Doos  method.  In  his  new  method,  Cressman  used  a  pre¬ 
liminary  field  which  was  usually  the  forecast  for  the  time  of  the 
observation  data.  Weighted  corrections,  based  on  the  new  observations, 
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were  then  applied  to  the  preliminary  field.  The  corrections  were 
defined  as  a  function  of  the  distance  (d)  from  grid  point  to  station. 

The  weight  factor  (W)  for  each  correction  is  given  by 

„  n2  -  d2  . 

W  =  n2  +  d2  tor  d  -  n 

=  0  for  d  >  n 

where  n  is  a  multiple  of  N,  the  grid  interval.  The  correction  procedure 
was  then  repeated  for  decreasing  multiples  of  N.  Cressman's  new  method 
of  objective  analysis  resulted  in  better  numerical  predictions  than 
those  predictions  based  upon  his  previous  method  of  analysis . 

Improvements  of  existing  techniques  have  been  developed,  most  of 
which  are  based  on  the  inclusion  of  additional  information  such  as 
surface  data,  vorticity,  geostrophic  approximation,  etc.  (see  Doos  and 
Eaton,  1957;  Johnson,  1957;  Sasaki,  1958;  Aubert,  1959;  Masuda  and 
Arakawa,  1962;  and  Teweles  and  Snidero,  1962).. 

Although  the  above  objective  analyses  give  satisfactory  results 
for  regions  of  sufficiently  dense  observation  stations,  it  remains  to 
develop  a  reliable  technique  which  will  operate  satisfactorily  over 
regions  of  sparse  data.  An  attempt  to  develop  such  a  technique  is 
presented  in  the  following  sections . 


2.  Basic  Spline  Theory 


The  problem  of  passing  a  snooth  curve  through  a  given  set  of 
points  (H)  has  been  solved  Mechanically  by  using  a  thin  elastic  strip 
to  define  the  curve  (Figure  la),  this  strip  is  coanmly  called  a  spline. 
Although  a  polynomial  ot  degree  H-l  could  be  determined  which  uould  also 
pass  through  the  same  set  of  points,  the  curve  defined  by  the  spline 
will  be  smoother.  Froa  the  theory  of  elasticity  it  can  be  shown  that 
a  spline  will  have  the  B'nimua  possible  strain  energy  (Love,  1934; 
Holladay,  1957),  that  is 

/K24s  =  Minima 

where  K  is  the  curvature  and  ds  is  the  arc  length.  Since  strain  energy 
is  a  Measure  of  the  smoothness  of  a  curve,  the  spline  Must  define  the 
smoothest  curve  for  a  given  set  of  points. 

If,  in  addition  to  knowing  the  points  through  which  the  spline 
passes,  the  slopes  are  prescribed  at  these  points,  the  spline  will 
define  a  new  curve  (Figure  lfcj  having  a  minimum  of  strain  energy  for 
the  new  set  of  constraints. 

The  concept  of  obtaining  a  smooth  curve  by  p.issing  a  spline  through 
a  given  set  of  points  nay  be  applied  to  the  analysis  cf  data.  For  this 
purpose,  the  spline  curve  between  any  two  data  points  will  have  to  be 
approximated  mathematically  by  a  polynomial  representation.  These 
piecewise  polynomials  will  then  have  to  be  joined  together  under  certain 
specified  constraints  like  tte  continuity  of  the  function  and  its 
derivatives  at  the  lata  points.  These  constraints,  of  course,  depend  on 
the  nature  of  the  phenomenon  under  investigation.  Fowler  and  Wilson 
(1966)  have  developed  a  method  of  determining  a  series  of  cubic  equa¬ 
tions  which,  when  spliced  together  with  continuous  slope  and  curvature 
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Fig-  la:  Spline  passing  through  points  A,  B,  C,  and  D. 


Fig.  lb:  Spline  passing  through  points  A,  B,  C,  and  D  where 
slopes  have  been  prescribed  at  the  points. 


c 

Fig.  1c:  Spline  curve  which  results  when  one  of  the  constraints 


(point  C)  has  been  removed. 
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at  the  junction  points,  approximate  a  spline  for  a  given  set  of  data 
and  slopes.  Pearce  (1968)  used  piecewise  quadratic  polynomials  1:0  allow 
for  the  sharp  but  real  changes  which  occur  in  representing  a  vertical 
wind  profile. 

Normally  it  is  not  desirable  to  fit  a  curve  exactly  through  all  the 
data  points,  and  so,  some  type  of  smoothing  is  performed  on  the  data. 

It  is  interesting  to  note  now,  that  if  one  of  the  constraints  in  Figure 
lb  is  removed,  the  resulting  position  of  the  spline  (Figure  lc)  is  one 
of  even  less  strain  energy.  That  is,  the  spline  defines  a  rmoother 
curve.  Thus  an  equation  which  approximates  the  new  spline  should  give 
a  smoother  representation  of  the  data.  However,  since  it  is  not  desir¬ 
able  to  completely  eliminate  any  data  points,  the  new  spline  could 
instead  be  used  to  determine  the  magnitude  and  direction  of  the  movement 
of  a  data  point  such  that  "controlled"  smoothing  may  be  done.  Fowler 
and  Wilson  (ibid.)  made  use  of  this  idea  to  smooth  their  series  of 
continuous  cubic  equations. 


# 


3.  Mathematical  Approximation  to  Splines 


The  mathematical  development  of  the  more  widely  used  spline  approx¬ 
imation,  the  cubic  spline,  is  described  below.  A  cubic  spline  fitting 
requires  a  general  third  degree  polynomial  between  any  two  points  of  a 
given  set  of  data.  Hence,  ten  initial  conditions  must  be  known.  If  the 
coordinate  system  is  translated  and  rotated  so  that  the  first  point  is 
at  the  origin  and  the  second  is  on  the  x'-axis,  the  general  equation  for 
a  third  degree  polynomial  reduces  to 

y'  -  Ax'3  +  Bx'2  +  Cx  +  D  (1) 

where  the  primes  indicate  the  new  coordinate  axes .  The  solution  to  this 
new  system  requires  only  four  initial  conditions.  Assuming  that  the 
coordinates  of  the  two  points  (endpoints)  are  known,  it  remains  to 
determine  two  additional  conditions.  These  are  obtained  in  the  follow¬ 
ing  manner.  Consider  Figure  2.  A  circle  is  fitted  to  the  set  of  points 
defined  by  the  endpoint  B,  and  the  two  adjacent  points,  A  and  C.  The 
derivative  of  the  resulting  equation  may  then  be  solved  for  the  slope 
at  B.  Similarly,  the  slope  at  E  is  determined  using  points  D,  E  and  F. 
The  four  initial  conditions  are  then: 


at  the  first  endpoint  x*  «  y'  ■  0  (2) 

slope  *  SI'  (3) 

at  the  second  endpoint  x'  »  d,  y*  ■  0  (4) 

slope  -  S2'  (5) 


where  d  is  the  distance  between  the  endpoints. 

Ihe  equation  for  the  slope  is  given  by  the  first  derivative  of  (1) . 

-  3Ax'  +  2Bx'  +  C  (6) 
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Fig.  2:  Determination  of  slopes  at  data  points  by  fitting  circles 
to  the  data. 


Fig.  3:  Normal  position  assumed  by  the  cubic  spline  when  fitted 


to  the  endpoints  of  a  three  data  interval. 
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The  system  consisting  of  equations  (1)  through  (6)  may  now  be  solved  in 
terms  of  SI',  S2'  and  d.  The  resulting  expression  for  y',  then,  is 

y.  ■  <&:.  +  S2-)  „,3  .  (2Siy.  S.2'1  +  sl.x.  (7) 


Due  to  the  rotation  of  the  coordinate  system,  the  endpoint  slopes  in  the 
original  system,  SI  and  S 2 ,  have  been  transformed  to  new  values,  SI'  and 
S2',  in  the  primed  system.  This  transformation  is  given  by: 


SI' 


SI  -  TR 
1  +  SI  •  TR 


S2 ' 


S2  -  TR 
1  +  S2  •  TR 


where  TR  is  the  tangent  of  the  angle  of  rotation  (0) . 

Now  consider  a  data  interval  to  be  the  interval  defined  by  the 
endpoints  of  three  successive  data  points.  A  cubic  equation  fitted  to 
the  endpoints  of  each  successive  interval  for  a  given  set  of  data  will 
usually  pass  between  the  curve  to  be  approximated  and  the  middle  data 
point  of  the  interval  (see  Figure  3).  This  suggests  an  iterative 
smoothing  procedure  wh1'  ch  will  result  in  the  convergence  of  the  data 
towards  the  curve . 

Starting  with  the  first  interval,  a  cubic  is  constructed  for  that 
interval.  The  middle  data  point  is  then  adjusted  towards  the  cubic. 
Since  the  location  of  the  curve  to  be  approximated  is  not  known  in  most 
cases,  the  middle  point  is  only  moved  some  fraction  of  the  total  dis¬ 
tance  between  the  point  and  the  cubic  curve.  The  procedure  is  repeated 
in  the  succeeding  intervals.  Figure  4a  demonstrates  this  procedure 
where  the  adjusted  points  are  indicated  by  the  higher  subscripts.  For 


cubic  from  A  to  C 

o  o  '■'0 

cubic  from  B  to  D 
o  o 

cubic  from  C  to  E 
o  o 


curve  to  be  approximated 
Fig.  4a:  Set  of  three  cubic  splines  approximating  a  curve, 


Fig.  4b:  Second  iteration  in  the  cubic  spline  convergence 
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example,  a  cubic  is  fitted  to  Aq  and  Co*  The  middle  point,  Bq,  is  then 
moved  toward  the  cubic  to  its  adjusted  position,  .  In  the  next  inter¬ 
val,  the  cubic  is  fitted  to  B0  and  Dq  which  adjusts  Cq  to  Cj .  Similarly 
Do  is  adjusted  to  Dj  in  the  interval  (Cq ,  Eq) .  Upon  completion  of  the 
data  adjustment  for  all  intervals,  the  entire  process  is  repeated  until 
satisfactory  convergence  is  achieved  or  the  benefit  of  additional 
smoothings  becomes  impractical.  Figures  4b  and  4c  depict  the  second 
and  third  iterations  in  a  series  of  three  smoothing  iterations  for  a 
set  of  five  data  points.  In  general,  spline  interpolation  converges 
rapidly  as  shown  by  Ahlberg  and  Nilson  (1963) ,  and  has  the  additional 
property  of  being  insensitive  to  round-off  errors. 


.  Two-Dimensional  Considerations 

The  results  of  approximating  curves  by  splines  have  been  so 
successful  (Curtis  and  Powell,  1966;  Walsh,  Ahlberg  and  Nilson,  19£2) 
that  it  would  seem  desirable  to  extend  the  theory  to  approximate 
surfaces.  Birkhoff  (1966)  developed  a  method  using  bicubic  spline 
interpolation  which  approximates  surfaces  on  a  rectangular  field  given 
the  data  value  (uj  at  all  grid  points,  the  normal  derivative  (3u/3n)  at 
the  boundary  grid  points  of  each  elemental  rectangle,  and  the  cross 
derivatives  (32u/3x3y)  at  the  four  corners  of  the  field.  In  practice, 
however,  most  data  fields  are  not  known  at  a  regular  grid  network  and 
the  problem  is  therefore  to  interpolate  the  known  data  to  the  grid 
points  and  smooth  the  resulting  fields. 

Since  circles  were  used  successfully  by  Fowler  and  Wilson  (ibid.) 
to  obtain  an  approximation  to  the  slope  at  each  point  in  their  curve 
fitting  routine,  it  would  seem  logical  to  fit  spheres  to  the  data  to 
approximate  the  surface  and  slope  in  a  surface  fitting  routine.  Since 
the  value  of  the  surface  is  desired  only  at  the  grid  points,  the  problem 
may  be  reduced  to  one  of  splicing  the  surfaces  of  the  spheres  together 
at  a  discrete  number  of  points. 

Given  an  arbitrarily  located  set  of  two  dimensional  data  points, 
the  points  are  ordered  according  to  increasing  values  of  the  x- 
coordinate.  A  grid  network  is  defined  (grid  length  is  arbitrary)  to 
cover  the  data  field  such  that  all  data  points  fall  within  the  network. 
Since  for  a  given  grid  point,  not  all  the  data  influence  the  value  at 
that  point,  what  shall  be  referred  to  as  the  "band  of  influence"  is 
defined  for  each  grid  line  (y  «  constant) .  The  vertical  plane  which  is 
coincident  with  the  grid  line  is  called  the  grid  plane. 
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For  a  given  grid  line  then,  those  data  points  which  fall  into  the 
band  of  influence  are  used  to  define  the  spheres.  The  surface  value  at 
each  grid  point  on  the  grid  line  can  then  be  determined  by  substituting 
the  coordinate  of  the  grid  point  into  the  equation  for  the  appropriate 
sphere.  The  equation  for  a  sphere  is  given  by: 

(x  -  a)2  +  (y  -  b)2  +  (z  -  c)2  =  r2  (8) 


Given  four  data  points,  (8)  may  be  solved  for  a,  t,  and  c,  the  coordi¬ 
nates  of  the  center  of  the  sphere,  and  r,  the  radius  of  the  sphere. 

For  a  given  grid  line,  y  =  y  ^  ,  which  passes  through  the  sphere,  and  grid 
point,  x  *  in  the  domain  of  the  sphere,  (8)  may  be  solved  for  the 
surface,  z^,  at  the  grid  point. 

In  order  to  define  the  surface  from  the  beginning  of  the  grid  to 
the  first  data  point,  the  first  three  data  points  in  the  band  of  influ¬ 
ence  are  used  to  quadratically  extrapolate  backwards  past  the  beginning 
of  the  grid.  Based  on  this  extrapolation,  three  dummy  data  points  are 
created.  The  requirement  on  the  location  of  these  dummy  points  is  that 
the  first  two  are  located  outside  the  grid  boundary  and  the  third  loca¬ 
ted  between  the  boundary  and  the  first  real  data  point.  The  reason  for 
this  will  be  explained  below.  Three  dunsny  points  are  also  determined 
at  the  other  end  of  the  grid  to  generate  data  from  the  last  real  data 
point  to  the  end  of  the  grid. 

The  surface  value  at  the  first  grid  point  along  a  grid  line  is 
determined  by  fitting  a  sphere  to  the  first  four  data  poincs  and  sub¬ 
stituting  the  grid  point  into  the  equation  for  the  sphere.  The  first 
grid  point  is  always  made  to  fall  between  the  second  and  third  data 
points.  If  the  succeeding  grid  point  also  falls  between  the  second  and 
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third  data  points,  its  surface  value  is  obtained  using  the  same  sphere. 
As  soon  as  the  succeeding  grid  point  falls  past  the  third  data  point,  a 
new  sphere  is  used  to  obtain  the  surface  value  for  the  new  grid  point. 
The  new  sphere  is  determined  from  the  last  three  data  points  of  the 
previous  sphere  and  the  next  data  point  in  succession  (see  Figure  5) . 
Thus,  surface  values  at  all  grid  points  represent  the  "midsection",  so 
to  speak,  of  the  data  which  defines  each  sphere. 

This  procedure  is  repeated  for  each  grid  line  until  an  approxima¬ 
tion  to  the  surface  has  been  determined  at  all  grid  points. 

The  resulting  data  values  at  the  grid  points  are  checked  for  error 
by  comparing  the  grid  point  value  to  the  mean  of  the  data  values  in  the 
band  of  influence  being  considered.  If  the  difference  between  a  grid 
point  value  and  the  mean  data  value  is  greater  than  two  standard  devia¬ 
tions  of  the  real  data,  the  grid  point  value  is  replaced  by  a  linear 
interpolation  between  the  two  nearest  acceptable  grid  values  along  the 
grid  line  (one  acceptable  grid  value  is  taken  from  each  side  of  the  re¬ 
jected  grid  value).  If  the  first  grid  point  is  rejected,  it  is  replaced 
by  the  mean  of  the  data  in  the  band  of  influence.  This  procedure  is  to 
some  extent  arbicrary,  since  the  mean  of  the  entire  data  field  may  be 
used  instead  of  the  mean  for  the  band  of  influence.  Also,  in  some 
analyses  it  may  be  desirable  to  restrict  the  magnitude  of  the  extreme 
values,  thus,  the  number  of  standard  deviations  of  the  real  data  is  left 
as  a  variable  input  parameter.  Since  it  is  possible  for  two  successive 
grid  point  values  to  occur  at  opposite  extremes  within  the  prescribed 
number  of  standard  deviations  of  real  data,  the  difference  between  two 
successive  grid  values  may  be  as  large  as  twice  the  number  of  standard 
deviations  wnich  was  prescribed.  Such  a  change  may  or  may  not  be 
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realistic,  sa  all  grid  poise  values  (except  those  it  *  *  *,)  are  sub¬ 
jected  to  an  additional  erro'.  check  which  is  based  upon  the  wiws 
persi8«iM«  slope  (KPS)  between  grid  points.  The  KPS  is  a  function  of 
the  variable  being  analyzed  and  is  therefore  an  inpat  parameter. 
Beginning  with  the  second  grid  point,  the  points,  x_,  along  the 
boundary,  y]t  of  the  grid  are  compared  to  the  grid  values  at  the  last 
acceptable  grid  point,  (x.,y,),  and  (x4 .y^).  If  the  slope  fit*  the 
point  (x.,y.)  to  any  two  of  the  neighboring  points  is  greater  than  KPS 
the  grid  point  value  at  (x^y^  is  replaced  in  the  sane  aanaer  as 

described  above.  For  the  points  (x.,y  ;  i  >  2,  1  <  j  <  j  >  each  grid 

1  J  rax 

value  along  the  grid  lines,  y ^ ,  is  compared  to  the  two  neighboring  grid 

values  (r  ,v.)  and  (x  ,y.  t>.  If  the  slope  iron  the  point  (x.,t.)  to 
ill  i  J~.  l  "j 

either  of  the  neighboring  points  is  greater  than  KPS,  the  grid  point 

value  is  replaced  as  described  above.  The  error  check  at  the  boundary, 

y  ,  is  analogous  to  the  error  check  at  the  boundary,  y  .  Errors 
Jmax  1 

occurring  near  the  end  of  a  grid  line  nay  not  have  an  acceptable  adja¬ 
cent  grid  value  in  the  positive  x-directien.  In  this  situation  a  second 
order  Taylor's  extrapolation  is  used  froa  the  last  three  acceptable 
grid  values  to  the  end  of  the  grid. 

Starting  with  the  first  grid  line,  y^,  each  grid  lice  is  snoot hed 
using  the  cubic  spline  routine  as  outlined  in  Section  3.  Upon  conple- 
tion  of  smoothing  in  the  x-direction,  the  entire  field  is  saootned  in 
the  same  manner  in  the  y-direction.  Repeated  smoothings  in  the  x  and  y 
directions  eliminate  directional  bias,  and  after  three  or  four  of  these 
bi-directional  smoothings,  the  grid  point  values  converge  to  a  surface. 


5.  Results 


In  order  to  properly  test  the  spline  Method  it  would  be  advanta¬ 
geous  to  Saw  an  exact  solution  for  the  surface.  Thus ,  in  addition  to  a 
qualitative  evaluation  of  the  results ,  a  quantitative  Measure  of  the 
error  could  he  determined.  This  would  have  the  additional  advantage  of 
being  able  to  compare,  quantitatively,  the  error  of  the  spline  objective 
analysis  to  the  error  of  other  methods. 

based  on  these  considerations,  a  test  was  constructed  using  a  solu¬ 
tion  for  the  stream  function  for  the  upper  level  of  a  two  layer  baro- 
clinic  spectral  general  circulation  model.  The  solution  was  for  the 
northern  hemisphere  and  could  be  interpreted  as  the  height  of  the  250  oh 
pressure  surface.  A  to *a 1  of  356  data  points  were  selected,  most  of 
which  are  established  radiosonde  stations  (see  Figure  6).  An  18  x  72 
grid  (grid  length  *  5*  latitude)  was  used  and  the  height  at  all  grid 
points  was  determined  iron  the  exact  solution.  The  exact  solution  is 
shown  in  Figure  7. 

The  results  of  applying  the  spline  technique  to  the  data  at  the 
radiosonde  locations  are  shown  in  Figure  S.  Another  objective  analysis 
technique  (Cressman,  1959)  was  applied  to  the  same  data  and  the  results 
are  shown  in  Figure  9.  The  standard  error  was  computed  for  both  tecb- 
niqves  with  the  result  that  the  error  for  Cressman’ s  technique  was  on 
the  order  of  twice  the  error  of  the  spline  technique.  The  percent  error 
for  the  splice  so lutioa  was  calculated  at  each  grid  point  and  is  shown 
in  Figure  6.  The  only  region  there  greater  than  101  error  occurs  is 
located  over  th~t  section  of  the  Pacific  ocean  where  there  is  at  xtreme 
lack  of  data  (see  Figure  6) -  The  nean  error  for  the  entire  northern 
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Fig.  6:  Data  point  locations  (marked  by  X)  and  percent  error  for 
the  spline  approximation  to  the  exact  solution.  Areas  not  shaded  are 
less  than  5Z  error. 


latitude 


LONGITUDE— 

Fig.  7:  Exact  solution  for  the  stream  function  for  the  upper  level 


of  a  two  layer  baroclinic  spectral  general  circulation  model. 
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Fig.  8:  Spline  approximation  to  the  exact  solution.  Standard 
error  =  147. 
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Fig.  9:  .  Cressman's  approximation  to  the  exact  solution.  Standard 
error  *  321. 
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hemisphere  was  approximately  3%  and  in  regions  of  good  data  coverage 
this  figure  usually  dropped  to  0  to  2 %. 

In  practical  weather  analysis,  quite  frequently  it  becomes  neces¬ 
sary  to  analyze  pressure  gradient  and/or  temperature  discontinuities 
(fronts).  Since  one  of  the  functions  of  most  objective  analyses  is  to 
smooth  the  data,  often  the  discontinuity  is  lost  completely  or  is 
smoothed  to  such  an  extent  that  its  frontal  characteristics  are  no 
longer  significant.  Since  it  is  usually  desirable  to  maintain  frontal 
characteristics,  a  test  was  constructed  to  determine  how  the  cubic 
spline  technique  performs  on  discontinuities.  The  temperature  data  for 
the  United  States  on  December  28,  1966  was  selected  for  analysis.  On 
this  date  a  particularly  sharp  cold  front  was  located  in  the  mid-section 
of  the  country.  Figures  9,  11,  and  13,  which  are  the  spline  analyses  of 
the  front  at  various  pressure  surfaces,  indicate  that  the  loss  of 
frontal  characteristics  appears  to  be  minimal,  even  at  the  850  mb  level 
where  the  front  is  most  intense.  Figures  10,  12,  and  14  are  the  corre¬ 
sponding  subjective  analyses.  A  comparison  of  objective  and  subjective 
analyses  shows  good  agreement  both  along  and  across  the  front.  It  is 
possible  that  the  frontal  characteristics  may  be  even  better  resolved  if 
a  quadratic  spline  is  used  as,  for  example,  by  Pearce. 

For  the  same  time  and  location  as  the  temperature  data,  the  wind 
and  height  data  were  also  analyzed  objectively  (Figures  15,  17,  19,  21, 
23,  and  25)  and  subjectively  (Figures  16,  18,  20,  22,  24,  and  26). 
Examination  of  these  and  the  previous  results  indicates  that  unless  a 
particular  feature  of  a  data  field  (for  example,  a  steep  gradient,  a 
maximum,  or  a  minimum)  is  "supported"  by  more  than  one  or  two  pieces  of 
data,  thu  feature  is  smoothed  rather  markedly.  More  explicitly,  the 
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wind  maximum  located  at  x  =  17 ,  y  =  16  on  Figure  16  is  defined  by  only 
one  piece  of  data.  On  the  corresponding  spline  analysis,  Figure  15, 
this  maximum  is  not  well  defined.  Likewise,  the  wind  minimum  at  x  =  45, 
y  =  12  on  Figure  20  does  not  show  up  well  on  the  corresponding  spline 
analysis,  figure  19.  Based  upon  one  piece  of  data,  the  500  mb  tempera¬ 
ture  minimtim  is  split  into  two  sections  in  Figure  12,  whereas  the  spline 
analysis  of  the  same  data  produces  only  one  minimum  region.  Of  course, 
in  most  situations  there  are  more  than  one  or  two  pieces  of  data  to 
define  the  significant  features  and  so  there  is  generally  good  agreement 
between  the  subjective  analyses  and  the  spline  analyses. 
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tig.  9:  Spline  analysis  of  the  850  mb 
temperatures  for  central  United  States,  00  GMT, 
December  28,  1966. 


Fig.  10:  Subjective  analysis  of  the  850  mb 
temperatures  for  central  United  States,  00  GMT, 
December  28,  1966.  Radiosonde  locations  are 
indicated  by  X. 
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Fig.  12:  Subjective  analysis  of  the  500  mb 
temperatures  for  central  United  States,  00  CMT, 
December  28,  1966.  Radiosonde  locations  are 
indicated  by  X. 


Fig.  13:  Spline  analysis  of  the  300  mb 
temperatures  for  central  United  States,  00  GMT, 
December  28,  1966. 


Fig.  14:  Subjective  analysis  of  the  300  mb 
temperatures  for  central  United  States,  00  CMT, 
December  28,  1966.  Radiosonde  locations  are 
Indicated  by  X. 


Fig.  16:  Subjective  analysis  of  the  850  mb 
vinda  for  central  United  States,  00  GMT, 

December  28,  1966.  Radiosonde  locations  are 
indicated  by  X. 


Fig,  18:  Subjective  analysis  of  the  500  mb 
winds  for  central  United  States,  00  GMT, 

December  28,  1966.  Radiosonde  locations  are 
indicated  by  X. 
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Fig.  19:  Spline  analysis  of  the  300  mb  winds 
for  central  United  States,  00  GMT,  December  28,  1966. 


Fig.  20:  Subjective  analysis  of  the  300  mb 
winds  for  central  United  States,  00  GMT, 

December  28,  1966.  Radiosonde  locations  are 
indicated  by  X. 


Fig.  2i:  Spline  analysis  cf  the  8 50  mb 
geopotential  heights  for  central  United  States, 
00  GKT,  December  28,  1966. 


Fig.  22:  Subjective  analysis  of  the  850  mb 
geopotential  heights  for  central  United  States, 
00  GMT,  Decenber  28,  1966.  Radiosonde  locations 
are  indicated  by  X. 


Fi*.  24:  Subjective  analysis  of  the  500  nb 
geopotential  heights  for  centra)  Uni ten  States. 
00  OfT,  Decaber  28,  196*.  Radiosonde  locations 
ate  indicated  by  X. 


Fij-  26=  Subjective  ualysis  of  the  300  mb 
K*«Fote«tiaI  bei^its  for  central  united  States, 
00  Off,  Breaker  2S,  1966.  latfioseade  locations 
are  indicated  br  I. 


6.  Conclusions 


Objective  analysis  by  the  spline  technique  appears  to  be  a  satis¬ 
factory  method  for  two  dimensional  data  analysis.  Analysis  of  regions 
with  poor  data  coverage  also  appears  to  give  satisfactory  results  except 
in  those  situations  where  the  features  being  analyzed  are  defined  by 
less  than  three  pieces  of  data.  The  magnitude  of  the  gradient  of  the 
data  to  be  analyzed  does  not  seem  to  have  any  undesirable  effects  on  the 
performance  of  the  technique,  provided  that  the  input  parameters  are 
properly  defined. 
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two-dimensional  data  field  is  described.  In  order  to  use  this  technique  it  is 
necessary  to  obtain  an  initial  approximation  to  the  data  at  the  grid  points.  This 
is  achieved  by  fitting  spherical  surfaces  to  the  data.  Bi-directional  Spline 
interpolation  was  then  applied  repeatedly  on  the  grid  point  estimates  of  the  data 
to  produce  convergence  to  the  true  surface. 

The  Spline  interpolation  technique  and  another  objective  analysis  technique 
developed  by  Cilchrist  and  Cressman  are  tested  against  an  exact  solution  and  the 
resulting  analyses  are  compared.  Real  temperature,  geopotential  height,  and  wind 
data  for  various  pressure  surfaces  are  analysed  by  the  Spline  method  and  the 
results  are  compared  to  the  subjective  analyses  of  the  same  data. 
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